Electron energy spectrum and density of states for 
non-symmetric heterostructures in an in-plane magnetic field 

A. Hernandez-Cabrer jj] and P. Accituno 
Dpto. Fisica Bdsica, Universidad de La Laguna, La Laguna, 38206-Tenerife, Spain 

F.T. VaskcQ 

Institute of Semiconductor Physics, NAS Ukraine, Pr. Nauki J^l, Kiev, 03028, Ukraine 

(Dated: February 6, 2008) 

Modifications of spin-splitting dispersion relations and density of states for electrons in non- 
symmetric heterostructures under in-plane magnetic field are studied within the envelope function 
formalism. Spin-orbit interactions, caused by both a slow potential and the heterojunction potentials 
(which are described by the boundary conditions) are taken into account. The interplay between 
these contributions and the magnetic field contribution to the spin-splitting term in the Hamiltonian 
is essential when energy amount resulting from the Zeeman and spin-orbit coupling are of the 
same order. Such modifications of the energy spectra allow us to separate the spin-orbit splitting 
contributions due to a slow potential and due to the heterojunctions. Numerical estimates for 
selectively-doped heterojunction and quantum well with narrow-gap region of electron localization 
are performed. 

PACS numbers: 72.25.-b, 73.21.-b 

I. INTRODUCTION 

Spin-orbit splitting of the energy dispersion relations for electrons in nonsymmetric quantum heterostructures has 
been theoretically considered during the past decades (see Ref. Q] for a review) . In bulk materials spin-orbit interaction 
appears both due to a slow-variable potential (related to the lattice constant) 2 and due to cubio* and linear— spin- 
dependent contributions to the effective Hamiltonian. Turning to the two-dimensional (2D) case, we can reduce the 
cubic contribution to a linear one after the replacement of the squared momentum by the quantized value due to 
confinement. It is still more important the fact that we have to take into account an additional spin-orbit splitting of 
the energy spectrum due to the interaction with abrupt heterojunction potentials (see Ref. 6 and discussion in Refs. 
20) • Such contribution is of a radically different kind with respect to the listed above because contributions from 
both sides of a slow confinement potential compensate each other 8 and, therefore, the spin-splitting of 2D states can 
not be obtained without a short-range potential contribution. To the best of our knowledge, the relative contributions 
from the bulk-induced mechanisms 2 *^ and from the heterojunctions have not been clarified experimentally in spite 
of a set of existing theoretical calculations 1 . In this paper we have examined the effect of an in-plane magnetic field 
on the electron energy spectrum and density of states and we have found that magneto-induced modifications of these 
characteristics are essentially different from the 2D spin-orbit interaction and Zeeman splitting. In principle, this fact 
makes possible an experimental verification of the above discussed contributions of the spin-orbit interaction. 

The effect of in-plane magnetic field on the energy spectrum in non-symmetric heterostructures occurs due to the 
Zeeman term in the electron Hamiltonian. This fact has been found i n Ref. shortly after initial considerations of 
the spin-orbit splitting in non-symmetric heterostructures (see Refs. I il l 1(11 11) . A number of peculiarities in transport 
phenomena for 2D systems under an in-plane magnetic field were also discussectii*i2ii£ and more complicated cases, 




FIG. 1: Band diagrams for selectivly-doped heterojunction (a) and quantum well (&). Dashed lines show electron ground state 
energy levels and thin curves correspond to wave functions; d is the QW's width, AE C<V are the band offsets, and m is the 
effective mass in the narrow-gap region, while mb and m± are the barrier effective masses. 
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such as quasi-lD transport or spin Hall effect (see Refs. ini^ or—, respectively) were recently considered. All this 
papers only take into account the mix between the Zeeman contribution and the effective 2D spin-orbit interaction. 
Here we perform the calculations based on the three-band Kane model** with non-symmetric boundary conditions, 
which is valid for narrow-gap heterostructures, and with a slow potential described self-consistently (see typical band 
diagrams in Fig. 1). The analysis of magneto-induced modifications of 2D energy spectra and corresponding density 
of states is presented for typical parameters of InGaAs/InAlAs-selectively-doped hetero junctions, and InGaAs and 
InSb-based quantum wells (QWs). 

The paper is organized as follows. In the next section we discuss the conduction c-band eigenstate problem and 
describe the electronic states in non-symmetric narrow-gap heterostructures using the boundary conditions for the 
wave functions at the interfaces in the parabolic approximation. In Sec. Ill we solve the 2D-eigenstate problem using 
the averaged transverse field approach. Numerical self-consistent calculations of the magneto-induced (caused by a 
magnetic field) modifications of the energy spectra and the density of states are presented in Sec. IV. Conclusions 
are given in the last section. 



II. BASIC EQUATIONS 

We start from the formulation of the eigenstate problem for the electronic state of the conduction c-band localized 
in the selectively-doped heterojunction or in the non-symmetric QW. All the above listed spin-dependent contri- 
butions are taken into account in the analysis. We also present the density of states which is connected to the 
photoluminescence excitation (PLE) intensity. 



A. Three-band Kane model 



The electronic states in narrow-gap heterostructures are described by the three-band Kane matrix Hamiltonian 16 



e z + (v • 7r), 



7T = p- -A, 

c 



(1) 



where the kinetic momentum, 7r, contains the vector potential A = (Hz, 0,0), H||OY is an in-plane magnetic field 
and p = (p,p z ) is written in the p, z-representation through the 2D momentum p. Here we have also introduced 
the diagonal energy matrix e z whose elements determine the positions of the band extrema (the energy values e cz , 
Shzi and Eiz correspond to electron, heavy- and light-hole extrema) and the interband velocity matrix v. The 6x6 
Hermitian matrix v is determined by the following non-zero matrix elements 7 : 
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and V is the characteristic interband velocity for the Kane model. 

Neglecting the small contributions from other bands we suppose in (1) that the inverse heavy hole effective mass 

is equal to zero. For such a case valence u-band components of the wave functions, ipp Z , are expressed in terms of 

fl 2) 

the c-band spinor with components ip pz ' an d the effective Schrodinger equation for this spinor takes the form (see 
Ref. ^3 where the case H = have been considered) : 
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where n± = (n x ± iTr y ) are the circular components of the in- plane momenta, 
neglected here the strain effect due to lattice mismatch^ and used e vz = Eh z 
in Eq. (3) the z-dependent effective mass, m z 
according to the relations: 



(3) 

For the sake of simplicity we have 
= Ei z . It is convenient to introduce 
the characteristic spin velocity, V z , and the effective g- factor, g z , 
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In such definitions Eq. (3) may be rewritten as the spinor eigenvalue problem: 



<e cz - E + y - + p z - — p z 

Zm z Zm z 

-V z [a x p] 2 + ^ B Ha y ] * p2 = 0, (5) 

( 1 21 a 

where spinor ^pj, is determined by the components 4i pz , /jb = \g\H/ (m e c) is Bohr magneton, and a is the Pauli 
matrix. The energy values of c- and w-band extrema for the narrow-gap (A) and wide-gap (B) regions take the forms: 

£ C z = U z , e vz = -e g + U z , (A) (6) 
e cz = AE C + U Z , e vz = -e g + AE v + U z , (B) 

and the corresponding band diagrams for QW and selectively-doped heterojunction were shown in Fig.l. Here the 
energy is counted from the bottom of c-band, e g is the gap, AE C and AE V are the band offsets for c- and w-bands, 
correspondingly. Here the slow potential U z should be determined from a self-consistent procedure. Thus, after 
substitution of (6) in Eq. (5) we have formulated the eigenstate problem for the spinor \l/p Z . 

Since AE C > e g in narrow-gap heterostructures, then the weak underbarrier penetration of wave function takes 
place for the electronic states with E <§; e g (parabolic band approximation). Because of this, one can neglect the 
longitudinal motion for the underbarrier region, z < for the single heterojunction case shown in Fig. 1(a), where the 
solutions take the form: 

* * pz =oe KZ , z < 0. (7) 

Here Hk = \/2rribAE c and Redetermines the scale of the underbarrier penetration of wave functions. We have 
also used in Eq. (7) the continuity conditions for eigenfunctions at heterojunctions. Since V z in Eq. (4) 
is proportional to (de vz /dz), then the integration (5) over heterojunction produces additional contributions to the 
boundary condition for flows. Such contributions are proportional to the band offset at the heterojunction AE V . 
Eliminating the underbarrier contributions from such equation by the use of the explicit expression (7) we obtain the 
third kind boundary conditions^: 

(fe*p*)U=o - iPo^pz=o ~ *X[* x p] z * pz=0 = 0. (8) 

Here the momenta P — hum/nib characterizes the underbarrier penetration and the parameter x determines the 
spin-orbit coupling due to the abrupt potential of the heterojunction 

2m b f s AE V 
and the right-side part is written for the approximation nib — m 



B. Eigenstate problem 

In the parabolic approximation, we describe the c-band electronic states using the above introduced spinor ^f pz . 
Below we consider not very strong magnetic fields (huj c <C v F H/d, where v F is the Fermi velocity and lo c is the 
cyclotron frequency) when (p x — eHz/c) 2 in Eq. (5) is replaced by p 2 and the isotropic kinetic energy is given by 
E p = (p 2 +Py)/(2m) and includes the effective mass m. Considering the low-doped structure case (if e g s F ) we 
rewrite the Schrodinger equation (5) for the narrow-gap region in the form: 

p2z - + U z + W z ) = E^ pzi 

W z = -V z [&xp] z + ^-fM B Ha y (10) 
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where U z is the self-consistent potential, as given in Eq. (14) below, and the effective mass is ^-independent. The spin 
velocity V z and the z-dependent contribution to the <?-factor are proportional to the transverse electric field dU z /dz: 

v _ h du z ^ ~ ^ f 1 + JL!?E±\ (11) 

4mEg dz ' z 2m \ e g dz J 

For the selectively-doped heterojunction case we consider Eq. (10) for the z > region with the boundary conditions 
(8) and with ^p Z ^oo = 0. For the case of a QW of width d, one can eliminate the underbarrier contributions in 
analogy to Eqs. (7-9) and, in addition to Eq. (10) for the region \z\ < d/2, the following third kind boundary 
conditions should be used: 



(0«*p»)I*=±| T *^±* p ,±| + iX± [a- x p] z * p i±i = 0. (12) 

Here momenta P± = W2m±AE±(m/m±) determine the scale of the underbarrier penetration of wave functions; the 
different values of AE± and m± take into account the differences of band offsets and effective masses, as it is shown 
in Fig. 1(6). The parameters \± — AE v /2e g determine the spin-orbit coupling due to the abrupt potential of the 
heterojunction according to Eq. (9). When m + = m_ , which corresponds to the symmetric case, we will take 
P± = P and x± = Xi as we can see horn (9). 

Thus, both interface potentials, which determine the spin-dependent contributions to (12), and the intrawell field, 
which determines the spin velocity V z in (11), are responsible for the spin-splitting of energy spectra. Since the 
Zeeman spin-splitting term only appears in the Eq. (10), the mentioned magneto-induced modifications of electron 
states due to these two contributions are different. 

The self-consistent numerical procedure for the eigenstate problem (10) involves the potential U z , which is obtained 
from the Poisson equation in the following form: 

47re 2 [ z 

U z = / dz (z - z )[n D (z ) - n e (z )] (13) 

^ J —00 

Here n D (z) is the 3D concentration of donors and e is the dielectric permittivity that we have supposed as uniform 

across the heterostructure. The electron density distribution is introduced through the electron dispersion relations 
E ap according to: 

»»(*) = E / (2^)2 *pi + ' *pM £ - ~ E °p), ( 14 ) 

where ^p Z + denotes the Hermitian conjugate of and a = ±1 refers to the two possible spin orientations. The 
Heaviside function, 0(x), appears here for the zero-temperature case. The Fermi energy, e F , is expressed through the 
total electron density, n 2D , defined as n 2D = J dzn e (z). Thus, e F depends on H for the fixed concentration case. 
The density of states is given by the standard formula: 

^=E/t^%-^p)- (15) 

In order to analyze p e one needs to solve the above-formulated eigenstate problem and to perform the integrations in 
Eq. (15). The density of states is connected to PLE intensity for the case of near-edge transitions, I PLE - Since the 
interband matrix element v cv do not depend on the in-plane quantum numbers, one obtain a 7 ! 19 : 

Iple ~ E l e ' V ™| 2<5 ( £ >c - £ A t , - hu) ~ PHAlj, (16) 
A C A„ 

where e is the polarization vector, Au> — to — e g /h and e g is the gap energy, which is renormalized due to the 
confinement effect. 



III. ANALYTICAL CONSIDERATION 



Before the numerical consideration we perform a simplified calculation of the eigenstate problem using the uniform 
transverse field approximation U z ~ eF ± z which is valid for non-doped QWs under an external modulating field F ± ; 
for heavy-doped structures F ± implies an averaged self-consistent transverse field. 
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A. Average field approach 



A simplification of the eigenstate problem (10) appears due the to ^-independent spin-orbit perturbation when 
W z in Eq. (10) becomes 



W — v [a x p] z + w H a y 



(17) 



with the averaged across structure characteristic spin velocity v ~ \e\F ± H/ '(4me s ) and the Zeeman splitting w H — 
(g/2)fi B H, where the z-dependent correction to the g-factor is neglected. The fundamental solutions ^ p2 for the Eq. 

(10), with the spin-dependent contribution (17), can be factorized as products of functions 'J'^ and ipi k ^ . Here the 
spinors for a = ±1 are determined from the eigenstate problem: 



(e p + W) * 



M 
P 



F <]><"> 



(18) 



whose solutions are giveni* by: 
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(vp + + w H )/iw p 

(vp_ + w H )/iw p 
1 



(19) 



"Pi 



j (VP X + W H ) 2 + (Vp y ) 2 



where p± — p x ± ip y and the energy values s± p describe the mix between internal spin-orbit interaction and Zeeman 
spin splitting. 

The two z-dependent fundamental solutions <fi z k ° ) (labeled below by k — a, b) are determined from the equation 



2m 1 



{ka) 
Z 1 



(20) 



written for the narrow-gap region. 

The general spinor solution ty pz is expressed through these fundamental solutions according to 



* p2 = ^(A+tf'+v+B+tp*^) 
^•- 1 ' + B_^- 1 ') 



(21) 



where the coefficients A±,B± are determined from the boundary conditions (12) or from (8) and the requirements 

*p Z ^oo = 0. 

After the substitution of the solution (21) for ip^^ = 0, which corresponds to the boundary condition at z — > oo, 
into (8) and the multiplication of this system by ^ ( p ) on the left side we rewrite the boundary condition at z — as 
follows: 



V+A + +x% +1)+ 



-p H 
p- 

-p, 
p- 



* P *=o = 0, 
* P *=o = 0, 



(22) 



with V± = {p z + iP ) f z a±) | z= o ■ In order to calculate the proportional to \ contributions we use here the solutions 
of the spin-dependent eigenstate problem (19). Thus, Eq. (22) can be transformed with the use of the relations 
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(23) 



to a simple linear system for A±. The dispersion relation, E av , is determined from the zero determinant requirement. 
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FIG. 2: (Color on line) (a) Dispersion laws for different h values. Solid, dashed and dot-dashed curves corresponds to the 
dimensionless momenta p y /mv =0, 0.25, and 1, respectively, (b) Density of states vs dimensionless energy for magnetic fields 
h —0 (solid), 0.5 (dashed), 1 (dotted), 2.5 (dash-dotted), and 5 (dash-dot-dotted). 



A similar transformation of the boundary conditions (12) for QWs, after substitution on (21), permit us to rewrite 
the boundary condition at z = ±d/2 in the form: 



-p+ 
V- o 



-p+ 
V- o 



z=±# 



(24) 



Thus, we have obtained a linear system for A± and B± and, therefore, E ap is obtained from the solvability condition 
for this system. 



B. 2D model 



First, let us consider the case of a heterostructure without spin-orbit contributions from heterojunctions, % = 0, 
when the dispersion relation is given by Eq. (19). For the zero magnetic field case, w H = 0, the energy w p is replaced 
by vp and e± p in Eq. (19) is transformed to the isotropic dispersion relation e„ v = e p + u\v\p. For this case the 
density of states was considered in Ref. 0- If w H ^ 0, the dispersion relation (19) becomes anisotropic as it is shown 
in Fig. 2 (a). We have used in this figure the dimensionless magnetic field h — w H /mv 2 . In order to represent a 
general case, valid for any structure, we have also used dimensionless energy axis. 

For the case of the strong magnetic field, when w p is replaced by w H , the Zeeman splitting effect appears to be 
dominant in the dispersion relation: e± p ~ e p ± w H . The spin-orbit splitting is mainly manifested as a shift of the 
dispersion paraboloids towards higher (lower) p x values for a = +1(— 1), respectively. On the other hand, Zeeman 
splitting is shown as a displacement to higher (lower) energy values depending on a. The anisotropy of the dispersion 
relation e ap appears to be essential for the region around the cross-point , e c = w^/2mv 2 . 

After the shift p x — > p x — w H /rriv and the integration over p with the use of the 5- function, the density of states 
is transformed from Eq. (15) into the integral with respect to the cosine of the in-plane angle, £ = cos#: 
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Here £± = v a ^ ± 



1 + 2e/hw H and ^ = Q — <j/h. A straightforward integration for the above-cross-point 



region e > e c gives as the exact relation: p e = p 2D , while for the below-cross-point region e < e c we plot p E /p 2D versus 
dimensionless energy 2e/mv 2 for different fields h. Figure 2(6) shows this density of states for different dimensionlcss 
magnetic fields. Two points deserve special attention. First, the typical e -1 / 2 singularity around e = 0, which occurs 
for h = 0, tends to disappears with increasing field. Second, abrupt steps corresponding to the filling of each level 
are shifted with the field intensity indicating a delay in the filling caused by the Zeeman splitting. It should be noted 
that the lowest (a = —1) level seems to be overfilled for low magnetic fields, i.e., p E /p 2D > 1/2- Electrons equally 
distributes between levels for magnetic fields higher than h = 10. 



IV. NUMERICAL RESULTS 

Here we perform numerical calculations based on the simplified consideration outlined in Sec. IIS, as well as the 
self-consistent solution of the eigenstate problem of Sec. IL4. The correspondent level-splitting at the Fermi energy 
level in a selectively-doped heterojunction and the densities of states in QWs are described. 

A. Selectively-doped heterojunction 

The eigenfunction of Eq. (20) for z > is written through the Airy Ai-function , Ai[z/l± — (E — e ap )/e±] with 
l ± = {/h 2 /2m\e\F 1 _ and e±_ = {h/l±) 2 /2m, if the zero boundary condition at z — > oo is applied. After substitution of 
this function into the system (22), one obtains the dispersion equation in the form: 
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After calculating this determinant and introducing the dimensionless function K{x) 
Eq. (26) into 



0. (26) 
— Ai'{x)/Ai(x) one transforms 
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ej_w p 



(27) 



where v s = xs±l±/h is the spin velocity due to the interface contribution. The spin-dependent dispersion relation, 
E ap , is given by the roots of this equation. 

Let us consider the 2D approach, when E is around the ground energy level, e a , determined by the equation: 
P l±/h + K(—e /e±) = 0. Expanding K[(e ap — E)/e±] over E — e Q — e p one can transform Eq. (27) into a quadratic 
equation and the dispersion relation takes the form E ap ~ e Q + e p + aAE p , where the splitting AE p of the energy 
spectrum is given by 



AE n 




vp z + w H p x 
w 



+ 



V S W H Py 



(28) 



with the heterojunction-induced spin velocity, v s = v s /K(—e /e±), and the internal spin velocity, v, introduced in 
Eq. (17). In Fig. 3 we plot these dispersion relations for the same parameters used in Fig. 2(a). As in Fig. 2(a), we 
have used dimensionless energy and momenta to represent a general case. Comparing Fig. 2(a) with Fig. 3 we can 
see the effect of the interface contribution as an enhancement of the splitting between parabolas a = ±1. 

For high electron concentration the approach of Eq. (28) is no longer valid. In order to obtain E for this case we 
need to solve Eq. (10) together with Eq. (13) by means of a self-consistent approach. To do that we have used the 
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FIG. 3: (Color online) Energy dispersion relations calculated with equation (28) for v a — v for different magnetic fields h. Solid, 
dashed and dot-dashed curves corresponds to the dimensionless momenta p y /mv =0, 0.25, and 1, respectively. 
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FIG. 4: (Color online) Dispersion relations calculated numerically for different magnetic field H and n,2D = 10 12 cm 2 . Solid 
line: p y — 0. Dashed line: p y — 0.25 mv. Dot-dashed line: p y — mv. 



Transfer Matrix Method^. Fi gure 4 shows dispersion relations obtained in this way for Ino.8Alo.2As/Ino.52Gao.4sAs 
with the electron density n 2D — 10 12 cm~ 2 , which corresponds to a Fermi energy ef = 209 meV. Using the standard 
parameters2i one obtains x — 0.303 and v = 1.13 x 10 6 cm/s, so that dimensionless field h = 1 corresponds to 
H = 0.23 T. Thus, we can compare panel for h = 1 in Fig. 3 with panel for H = 0.23 T in Fig. 4. We can see 
that numerically calculated results show a similar behavior than approximation (28), although, for the material under 
consideration, splitting is qualitatively less and shift along momentum is larger in Fig. 4. 

The spin splitting of the dispersion relations at Fermi energy is directly connected to the Shubnikov-de Haas (SdH) 
oscillations. Since AE p is nearly isotrop ic over p-plane, we have also calculated AE = AE p \ p=PF for different 
magnetic fields, see Fig. 5. Following Refs. Il2ll3l the spin splitting is related to the modulation of the SdH oscillation 
amplitude according to cos(irmAE / h\e\H) . For the structure under consideration, experimental valued is AE = 11.4 
meV for H = 0.8 T and n 2D — 10 12 cm~ 2 , which corresponds to AE = 9.86 meV, calculated numerically for the 
same magnetic field and density. Figure 5 shows a very slight dependence of AE on the magnetic field for x = 0.303 
being stronger this dependence for the case x — 0- Thus, interface contributions should be detectable through the 
SdH oscillations amplitude. 



9 



10.0 




FIG. 5: The level splitting energy, AE vs magnetic field for \ = 0.303 and X = 0. Note that the upper AE is almost 
H-independent. 
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FIG. 6: (Color online) (a) Dispersion laws for InGaAs/InAlAs QW under different magnetic fields. p y — (solid line) and 
p y = mv (dashed line), (b) Density of states for H — T (solid line), 0.5 T (dashed line), 1 T (dot line), and 2.5 T (dot-dashed 
line). 



B. QW under field F_ 



Next, we consider the non-doped symmetric QW (with P± = P a and x± — X) subjected to a transversal field Fx, 

when the eigenstate functions of Eq. (20), ip[ a ' b \ are written through the Ai- and Bi-functions. The system (24) 
gives the dispersion relations equation as the determinant of a 4 x 4 matrix, which differs from Eq. (26) due to the 
replacement of V± and ^ by the 2x2 matrices: 



( \ (k,±l)\ 

[Pzfz J 



=±d/2 



_ ■ r, (fc,±l) 
T *Po<P ±d/2 ■ 



(acr) (bcr) 
(air) (6<t) 

f-d/2 V-d/2 



(29) 



Once again, after the expansion near the ground energy, e a , one can 



transform the determinant in a similar way to transformations of Eqs. (26, 27). Numerical solutions of these equations 
for InGaAs- and InSb-based QWs are represented in Figs. 6(a) and 7(a). We have considered 100 A wide InGaAs- and 
InSb-based QWs using the parameters of Refs. I2HI22T We have also applied a field of 120 kV/cm, which corresponds 
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FIG. 7: (Color online) (a) Dispersion relations for InSb/InAlSb QW under different magnetic fields. p y — (solid line) and 
p y — rav (dashed line), (b) Density of states for H = T (solid line), 0.25 T (dashed line), 0.5 T (dotted line), and 1 T 
(dot-dashed line). 



to spin velocity v = 1.04 x 10 6 cm/s, for the InGaAs QW, and a field of 100 kV/cm, spin velocity v = 8.2 x 10 6 cm/s, 
for the InSb. 

Comparing Figs. 6(a) and 7(a) with Fig. 4 we can see that, for similar spin velocities and magnetic fields, the 
splitting between levels is bigger whereas the shift in the p x direction is smaller for the QW case with respect to 
the selectively-doped structure. Thus, the effect of interface contribution seems to be more pronounced for the QW 
case. This effect is enhanced as the width of the wells diminishes. Due to the greater spin velocities of InGaAs and 
InSb-based QWs, caused by the different effective masses and energy gaps, it is necessary the use of lesser magnetic 
fields to get a similar effect. Both for QWs and selectively doped structures, interface contributions are opposite to 
the intrinsic spin-orbit coupling effect. 

Since the density of states is proportional to the PLE intensitjii^, as shown in Eq. (16), it is interesting to study 
the shape of p e , shown in Figs. 6(6) and 7(6). The effect of the interfaces is manifested in p £ as a delay in the 
quenching of the e~ ^-singularity corresponding to the zero-field case. Although the singularity no longer exists for 
H 7^ 0, a peak still remains at the energy value of the bands anticrossing. For the case which includes interface 
contributions, this peak gets wider and p e tends slower to p 2D for high energy limit. Thus, the PLE technique is of 
great interest to analyze the interface contributions 23 . It should be notice that a = ±1 levels becomes equally filled 
for magnetic fields beyond 5 T. 

V. CONCLUSION 

In this paper we have examined the electron states in narrow-gap non-symmetric heterostructures under an in- 
plane magnetic field. The cigenstate problem was formulated in the framework of the three-band Kane model with 
non-symmetric boundary conditions. We have found that the mechanisms of mixing between the Zeeman term in 
the Hamiltonian and the two kinds of spin-orbit coupling contributions (from a slow field and from heterojunctions) 
are essentially different. Numerical estimates for typical parameters of InGaAs/InAlAs and InSb/InAlSb structures 
demonstrate the essential magneto-induced modifications of the energy spectra under magnetic field strength of the 
order of Tesla. 

Let us discuss some possibilities for the experimental verification of the energy spectra modifications obtained here. 
The magnetotransport measurements of Shubnikov-de Haas oscillations under nearly in-plane magnetic fields^^^ 
(when a quasi-classic quantization of the dispersion relations for the transverse component of the magnetic field 
is possible) provide a direct information about in-plane magneto-induced modifications of electron energy spectra. 
The results of Sec. IV.A are in agreement with Ref. (3 but more measurements (for different in-plane fields and 
concentrations) are necessary in order to separate the intrinsic and junction-induced contributions. Another way of 
looking for these peculiarities of energy spectra is the mid-infrared PLE spectroscopjiS 3 . when interband transitions 
are modified under in-plane magnetic field. PLE intensity provides direct information of the energy spectra because 
it is directly connected to the density of states, as mentioned above. A comparison between SdH oscillations and PLE 
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measurements with a precision about 1 meV (it should be possible for a high-quality structure) provides more data 
on mechanisms of spin-orbit interaction in narrow-gap structures. 

Thus, we have shown the essential modifications of the electron energy spectra in non-symmetric narrow-gap QWs 
under in-plane magnetic fields. We suggest that measurements of the magneto- induced contributions to optical 
and transport properties would be an useful method for an experimental verification of the heterojunction-induced 
contributions to the spin-orbit interaction in the narrow-gap heterostructures under investigation. 
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